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ABSTRACT 



Context. We report on the analysis of XMM-Newton observations of RX J1347.5-1 145 (z=0.451), the most X-ray-luminous galaxy 
cluster. 

Aims. We present a detailed total and gas mass determination up to large distances (~ 1.7 Mpc), study the scaling properties of the 
cluster, and explore the role of AGN heating in the cluster cool core. 

Mefhods. By means of spatially resolved spectroscopy we derive density, temperature, entropy, and cooling time profiles of the intra- 
cluster medium. We compute the total mass profile of the cluster in the assumption of hydrostatic equilibrium. 
Results. If the disturbed south-east region of the cluster is excluded from the analysis, our results on shape, normalization, scaling 
properties of density, temperature, entropy, and cooling time profiles are fully consistent with those of relaxed, cool core clusters. 
We compare our total and gas mass estimates with previous X-ray, tensing, dynamical, and SZ studies. We find good agreement 
with other X-ray results, dynamical mass measurements, weak tensing masses and SZ results. We confirm a discrepancy of a factor 
~2 between strong tensing and X-ray mass determinations and find a gross mismatch between our total mass estimate and the mass 
reconstructed through the combination of both strong and weak tensing. We explore the effervescent heating scenario in the core of 
RX J1347.5-1 145 and find support to the picture that AGN outflows and heat conduction are able to quenching radiative cooling. 
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1. Introduction 

. X-ray observations of the diffuse Intra-Cluster Medium (ICM) 
' in clusters of galaxies are a particularly rich source of infor- 
mation for understanding the formation of large scale structure 
and the physics of clusters. As they are the last manifestation of 
hierarchical clustering, whose history depends strongly on cos- 
, mol ogy, galax y clusters are key objects for cosmological studies 
■ (see lVoil 12005 , for a review). Since the evolution of the ICM is 
' mainly driven by the gravity of the underlying dark matter halo, 
clusters are expected to show similar properties when rescaled 
with respect to their total mass and formation epoch. However, 
deviations from self-similarity are expected under the effect of 
more complex physical processes, beyond gravitational dynam- 
ics only, which affect the thermodynamical properties of the dif- 
fuse ICM (e.g. lEvrard & Henrvlll991l: iBrvan & Normanlll998t 
iBorgani etalj|2002 l. and references therein). It is therefore es- 



sential to investigate whether galaxy clusters obey the expected 
scaling relations, which are the foundation to use these virialized 
objects as cosmological probes. The first important step in this 
context is to find a proxy for an accurate determination of the 
cluster mass. 

The galaxy cluster RX J1347.5-1145 (z=0.451) is an ex- 
ceptional object in many aspects. It is the most X-ray-luminous 
cluster known to date (Lx = 6x10'*^ erg s"' in the [2-10] keV 
energy range) with a very peaked surface brightness profile and 
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hosts a strong cooling flow in its center with nominal ma ss ac- 
cretion rate of ~ 1900 Moyr ^ dOitti & Schindleill2004 . The 
cluster is dominated by two cD galaxies which are separated by 
about ~ 18" along the east- west direction, the X-ray emission 
being centered on the western one. Although this is unusual for 
strong cooling flow clusters, the optical spectrum of the western 
Brightest Cluster Galaxy (BCG) indicates that it hosts an active 
galactic nucleus (AGN), with typical emission lines of giant el- 
lipticals at the center of cooling flow clusters dCohen & KneibI 
2002). More striki ng is a recent discov ery made with Chandra 
dAUen et alj|2002bl) and XMM-Newton dGitti & Schindleiil2004l) 
of a region with hot, bright X-ray emission located at ~ 20 
arcsec from the central emission peak in south-east direction. 
Submillimeter observations also detected a very deep SZ decre- 
ment in the south-east regi on of the cluster (iKomatsu et alj 19991 : 
iPointecouteau et ani2001h . These results were interpreted as in- 
dications of a subcluster merger in an otherwise relaxed, massive 
cool core cluster, pointing to a complex dynamical evolution of 
the system. Furthermore, RX J1347.5-1 145 is a powerful gravi- 
tational lens and mass reconstructions based on weak and strong 
lensing analyses have been performed JSchindler et al. 1 995t 
iFischer & Tvsonlll997l:[Sahu et alJll998tlBradac et al.ll2005bl) . 

With a detailed study of the properties of the ICM in this 
cluster it is thus possible to address many key issues on both 
dynamical and non-gravitational processes in galaxy clusters. 
A great advantage of observing RX J1347.5-1145 with XMM- 
Newton is that important quantities derived for the undisturbed 
cluster (i.e., with the south-east quadrant excluded) such as the 
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azimuthally averaged ICM density and temperature profiles can 
be computed up to a large distance from the center (~ 1730 
kpc). The measurement of cluster temperature gradients at large 
distances is also crucial for determining the total gravitational 
masses and in turn the gas mass fraction of clusters. A pre- 
cise determination of the total mass at large radii allows an es- 
timate of the virial radius of the object without much e xtrapo- 
lation of the universal NFW dark matter profile (Navar ro et alj 
119961) . The virial radius can then be used to study the scalings of 
the temperature and entropy profiles and a fair comparison be- 
tween predictions of numerical simulations and observations can 
be performed. Currently, the two most promising techniques for 
obtaining accurate determinations of cluster masses are X-ray 
observations, by deprojection of X-ray surface brightness com- 
bined with spectroscopic determination of the cluster tempera- 
ture, and gravitational lensing, through either strong lensing fea- 
tures or statistical distortions of background objects (weak lens- 
ing). The mass estimates inferred with these two methods can be 
quite inconsiste nt, particularly in the case of strong lensing (e.g. 
IWu et al.|[T998 ^, and references therein). In contrast to the X-ray 
technique, the gravitational lensing method is essentially free of 
assumptions on the nature and the dynamical state of the grav- 
itating material. In partic ular, the X-ray m ethod can be aff'ected 
strongly during mergers (ISchindledl 19961) and in the inner clus- 
ter region where a strong interaction between the central AGN 
and the ICM is present (e.g., Birzan et al. 2004), as in these cases 
deviations from the assumptions of hydrostatic equilibrium and 
spherical symmetry are expected. Since both the total mass pro- 
file derived from X-rays and the total mass distribution derived 
from gravitational lensing are available for RX J1347.5-1 145, 
a comparison between them is possible thus providing impor- 
tant insights on this issue. Furthermore, the presence of gas with 
short cooling time in the cluster core offers the opportunity to 
explore gas heating processes such as AGN heating, which have 
become increasingly popular since the failure of standard cool- 
ing flows models. 

In this paper, by starting from the results of morphological 
(Sect. [3]l and spectral (Sect. |4| analyses of XMM-Newton ob- 
servations of RX J1347. 5-1145 (Sect.|2]i, we present a detailed 
study of the cluster mass distribution (Sect. |6]l, and discuss its 
comparison with the mass profile derived from previous stud- 
ies (Sect. |2ll. We also study the scaling properties of the cluster 
(Sect. |5] and |6]l and explore the role of AGN heating in the clus- 
ter cool core in the context of the eff'ervescent heating scenario 
(Sect. ID. RX J 1 347.5- 1145 (hereafter RX J 1 347) is at a redshift 
of 0.45 1 . With a Hubble constant of Hq - 70 km s"' Mpc"' , and 
= 1 - = 0.3, the luminosity distance is 2506 Mpc and 
the angular scale is 5.77 kpc per arcsec. 

2. Observation and data preparation 

RX J1347 was observed by XMM-Newton in July 2002 during 
rev. 484 with the MOS and pn detectors in Full Frame Mode 
with THIN filter, for an exposure time of 37.8 ks for MOS and 
33.2 ks for pn. We use the SASv6.0.0 processing tasks em- 
chain and epchain to generate calibrated event files from raw 
data. Throughout this analysis single pixel events for the pn data 
(PATTERN 0) are selected, while for the MOS data sets the 
RATTERNs 0-12 are used. The removal of bright pixels and hot 
columns is done in a conservative way applying the expression 
(FLAG==0). To reject the soft proton flares we accumulate the 
Ught curve in the [10-12] keV band for MOS and [12-14] keV 
band for pn, where the emission is dominated by the particle- 
induced background, and exclude all the intervals of exposure 



time having a count rate higher than a certain threshold value 
(the chosen threshold values are 0.15 cps for MOS and 0.22 cps 
for pn). The remaining exposure times after cleaning are 32.2 ks 
for MOSl, 32.5 ks for M0S2 and 27.9 ks for pn. Starting from 
the output of the SAS detection source task, we make a visual se- 
lection on a wide energy band MOS & pn image of point sources 
in the field of view (hereafter FOV). Events from these regions 
are excluded directly from each event list. 

The background estimates are obtained using a blank-sky 
observation cons isting of several h igh-latitude pointings with 
sources removed dLumb et"aLll2002l) . The blank-sky background 
events are selected using the same selection criteria (such as 
RATTERN, FLAG, etc.), intensity filter (for flare rejection) and 
point source removal used for the observation events; this yields 
final exposure times for the blank fields of 365 ks for MOSl, 
350 ks for MOS2 and 294 ks for pn. Since the cosmic ray in- 
duced background might slightly change with time, we compute 
the ratio of the total count rates in the high energy band ([10-12] 
keV for MOS and [12-14] keV for pn). The obtained normal- 
ization factors (0.992, 1.059, 1.273 for MOSl, MOS2 and pn, 
respectively) are then used to renormalize the blank field data. 
Furthermore, the blank-sky background files are recast in order 
to have the same sky coordinates as RX J 1347. The background 
subtraction (for spectra and surface brightness profiles) is per- 
formed as described in full detail in lArnaud et alj (l2002i) . This 
procedure consists of two steps. In a first step, for each prod- 
uct extracted from the observation event list, an equivalent prod- 
uct is extracted from the corresponding blank-field file and then 
subtracted from it. This allows us to remove the particle back- 
ground. However, if the background in the observation region is 
different from the average background in blank field data, this 
step could leave a residual background component. The residual 
background component is estimated by using blank field sub- 
tracted data in a region free of cluster emission and then sub- 
tracted in a second step from each MOS and pn product. 

The source and background events are corr ected for vi- 
gnettin g using the weighted method described in lArnaud et al.l 
('2001'), the weight coefficients being tabulated in the event list 
with the SAS task evigweight. This allows us to use the on-axis 
response matrices and effective areas. Unless otherwise stated, 
the reported errors are at 90% confidence level. 

3. Surface briglitness profile 

Previous Chandra and XMM-Newton observations of RX J 1347 
revealed the presence of a hot and bright X-ray subclump visible 
to the south-east (SE) of the main X-ray sur face brightness peak 
(lAUen et alJl2002bl: iGitti & Schindleiil2004 . On the other hand, 
the data excluding the SE quadrant (hereafter "undisturbed clus- 
ter") show a regular morphology, indicating a relaxed state. We 
are interested in determining the characteristic properties of the 
cluster in order to perform studies of mass profiles and scaling 
relations as it is usually done for relaxed clusters. The disturbed 
SE quadrant is thus masked in the following morphological anal- 
ysis. 

We compute a background-subtracted, vignetting-corrected, 
radial surface brightness profile in the [0.3-2] keV energy band 
for each camera separately. For the pn data, we generate a list 
of out-of-time event^ (hereafter OoT) to be treated as an ad- 
ditional background component. The effect of OoT in the cur- 

' Out-of-time events are caused by photons which arrive while the 
CCD is being read out, and are visible in an uncorrected image as a 
bright streak smeared out in RAWY. 
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Fig. 1. Background subtracted, azimuthally-averaged radial sur- 
face brightness profile in the [0.3-2] keV range of the data ex- 
cluding the SE quadrant (undisturbed cluster). The best fit /3- 
model fitted over the ~ 350 - 1730 kpc region is over-plotted as 
a dashed line (model SO in Table [TJ. When extrapolated to the 
center, this model shows a strong deficit as compared to the ob- 
served surface brightness. The solid line shows the best fit dou- 
ble yS-model fitted over the whole region (model DD in Table 



rent observing mode (Full Frame) is 6.3%. The OoT list is pro- 
cessed in a similar way as done for the pn observation event file. 
The profiles for the three detectors are then added into a sin- 
gle profile, binned such that at least a signal-to-noise ratio of 3 
is reached. The cluster emission is detected up to Ro^t - 1-73 
Mpc (~ 5 arcmin). The surface brightness profile of the undis- 
turbed cluster, shown in Fig.[T] is fitted in the CI AO tool Sherpa 
with various parametric models, which are convolved with the 
XMM point spread function (PSF). The overall PSF is obtained 
by adding the PSF of each camera (Ghizzardi 2001), estimated 
at an energy of 1.5 keV and weighted by the respective clus- 
ter count rate in the [0.3-2] keV en ergy band. A single y6-model 
dCavaliere & Fusco-Femianolll976h : 



^ r 



2 V -3/3+0.5 



(1) 



is not a good description of the entire profile (model SG in Table 
(TJ and a fit to the outer regions (350 kpc^ r^ 1730 kpc) shows a 
strong excess in the center as compared to the model (see Fig.lTJ. 
The centrally peaked emission is a strong indication of a cool- 
ing flow in this cluster We find that for 350 kpc ^ r < 1730 
kpc the data can be described by a y6-model with a core radius 

- 307 + 9 kpc and a slope parameter /3 = 0.86 ± 0.02 (3cr 
confidence level). The single yS-model functional form is a con- 
venient representation of the gas density profile in the outer re- 
gions, which is used as a tracer for the potential. The parame- 
ters of this best fit are thus used in the following to estimate the 
cluster gas and total mass profiles in the region where the single 
yS-model holds (see Sect.|6]l. 

We also consider a double isothermal yS-model in the form: 



1 + 



-3/3,+0.5 



(2) 



where / = 1,2, and find that it can account for the entire profile 
(see Fig. [Til. The best fit parameters are =39+1 kpc, jSi - 
0.62 ± 0.01, rc,2 = 386 ± 17 kpc,/32 = 1.01 + 0.05. By assuming 



a common /? value we find: rc i =241+7 kpc, rc_2 =47 + 2 kpc, 
P = 0.76 + 0.01 (see Table[B. 



4. Spectral analysis 

Throughout the analysis, a single spectrum is extracted for each 
region of interest and is then regrouped to reach a significance 
level of at least 25 counts in each bin. The d ata are modeled using 
the XSPEC code, version 11.3.0 ( Arnaudl l996l) . Unless other- 
wise stated,the relative normalizations of the MOS and pn spec- 
tra are left free when fitted simultaneously. We use the follow- 
ing response matrices: ml_439_im_pall_vl . 2 . rmf (MOSl), 
m2_439_im_pall_vl . 2 . rmf (MOS2), epn_f f 2(8_sY9 . rmf (pn). 

4.1. Global spectrum 

For each instrument, a global spectrum is extracted from all 
events lying within 5 arcmin to the cluster emission peak. We 
test in detail the consistency between the three camera by fit- 
ting separately these spectra with a mekal model (with the 
redshift fixed at z=0.451) absorbed by a column density in- 
cluded in the tbabs mo del (fixed at the nom inal galactic value 
A^H = 4.85 X 10^° cm ^ iDickev & Lockm anl990). Fitting the 
data from all instruments above 0.3 keV leads to inconsistent 
values for the temperature derived with the MOS and pn cam- 
eras: kT = 12.2+JJ-^ keV (MOSl), IQ.A^^l keV (MOS2), 9.3+°-^ 
keV (pn). We then perform a systematic study of the effect of 
imposing various high and low-energy cutoffs, for each instru- 
ment. Good agreement between the three cameras is found in 
the [0.8-10.0] keV energy range {kT 



1 1.2+°-^ keV for MOSl, 



10.0+°^ for MOS2, 10.2+JJ 4 for pn). We therefore perform the 
spectral analysis in this energy range. The combined MOS+pn 
global tempe rature, in keV, and metalli city, as a fraction of the 
solar value (1 Anders & Grevess^ Il989l) derived from the best 
fit OrVdof = 2717/1697) are respectively: kT = 10.4!°;3 keV, 
Z = 0.25;^QQ3 Zq. The unabsorbed luminosities in this model 
(estimated from the average of the fluxes measured by the three 
cameras after fixing A^h = 0) in the X-ray ([2.0-10.0] keV) and 
bolometric band are respectively: Lx - 6.2 ± 0.2 x 10^^ erg s"', 
^bol ~ 13.5+0.4x10^^ erg s"', where the errors are given as half 
the difference between the maximum and the minimum value. 



4.2. Spatially resolved spectra 

As done for the morphological analysis, for the spectral analy- 
sis we separate the SE quadrant containing the X-ray subclump 
from the rest of the cluster The data of the undisturbed cluster 
are divided into the following annular regions: 0-30", 30"- 1', 1'- 
1'.5, r.5-2', 2'-3', 3'-5'. The spectra are modeled using a sim- 
ple, single-temperature model (mekal plasma emission code in 
XSPEC) with the absorbing column density fixed at the nominal 
Galactic value. The free parameters in this model are the tem- 
perature kT, metallicity Z (measured relative to the solar values, 
with t he various elements assum ed to be present in their solar 
ratios, lAnders & Grevess"3 1 19891) and normalization (emission 
measure). The best-fitting parameter values and 90% confidence 
levels derived from the fits to the annular spectra are summarized 
in Table |2] 

4.3. Deprojectlon analysis 

Because of projection effects, the spectral properties at any point 
in the cluster are the emission-weighted superposition of radi- 
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Table 1. Results from fitting the surface brightness profile of the undisturbed cluster in different radial intervals [/?in-/?out]- The 
single and double jS-models used for the fitting are given by Eqs.[T]and|2l respectively. They are indicated with: SG (Single jS-model, 
fitted in the Global radial range), SO (Single y6-model, fitted in the Outer region), DD (Double j6-model, with Different /3 values), 
DE (Double j6-model, with Equal /3 values). The quoted errors are at 3cr confidence level. 



Model 


(arcmin) 


(kpc) 


(cts/s/arcmin^) 


P, 


(arcmin) 


(kpc) 




SG: single /? 
SO: single /? 
DD: double /3 
with j8i p-, 
DE: double /? 
with j8i = /32 


0.0-5.0 
1.0-5.0 
0.0-5.0 

0.0-5.0 


0-1731 
346-1731 
0-1731 

0-1731 


18 94+0-86 

18.12!j- 
0.96^- 


0.861!«|g 

6 16+™ 

1 010+""'^' 

Q7g, +0.012 
"•'"^-0.011 


0.8876!°;JJi^ 
0.1138tS 

1 1 1 4s: +0.0483 

1360+''™ 
6968+'' ™ 


52+j 
307~+^ 

40+1 
386+J2 

47!^ 

24 1+:^ 


1620/129 (12.56) 
109/87 (1.25) 
258/111 (2.32) 

289/112 (2.58) 



Table 2. Results of the spectral fitting in concentric annular 
regions in the [0.8-10.0] keV energy range obtained by fix- 
ing the absorbing column density to the Galactic value (A^h = 
4.85 X 10^'^cm"^). The temp erature (in keV) and metal licity (in 
fraction of the solar value, lAnders & Grevess^ Il989h are left 
as free parameters. The data of the SB quadrant are excluded 
(undisturbed cluster). 



Radius 


source counts 


kT 


Z 


;rVdof 


(kpc) 


(MOS+pn) 


(keV) 


(Z0) 




0-173 


46719 


q :!+0.J 
-0,3, 


31+0.05 


914/964 


173-346 


18377 


12 5+'-' 


16+''''' 

-flo; 


573/546 


346-519 


8733 


22+'' " 
-fi.15 


288/295 


519-692 


4331 


9 4+i-r 

9.8+j| 
7.3!« 


13+''''* 


201/178 


692-1038 
1038-1731 


4092 
2742 




315/229 
572/383 



ation originating at all points along the line of sight through 
the cluster. To correct for this effect, we perform a deprojection 
analysis by adopting the XSPEC pro jet model. Under the as- 
sumption of ellipsoidal (in our specific case, spherical) shells 
of emission, this model calculates the geometric weighting fac- 
tor according to which the emission is redistributed amongst the 
projected annuli. 

The deprojection analysis is performed by fitting simulta- 
neously the spectra of the three cameras. The results are re- 
ported in Table [3] We also calculate the electron density «e from 
the estimate of the Emission Integral EI - J n^n-pdV given by 
themekal normalization: 10"''*£'//(47r[DA(l +z)]^). We assume 
He = 1.2023«p in the ionized intra-cluster plasma. 

5. Radial profiles 

5.1. Temperature 

The deprojected temperature profile derived in Sect. l4.3l is shown 
in Fig. 121 where we also show the projected profile for compar- 
ison. As expected, the deprojected central temperature is lower 
than the projected one, since in the projected fits the spectrum 
of the central annulus is contaminated by hotter emission along 
the line of sight. We also not e that the projected temperature 
profile measured by Chandra d Allen et al] l2002lj) is systemat- 
ically slightly higher than that measured by XMM-Newton, al- 
though the general trend observed by the two satellites is consis- 
tent Tcitti & Schindler 2004). 

The temperature profile of RX J 1347 exhibits the shape 
characteristic for cool core clusters: the temperature declines 
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250 500 750 1000 1250 1500 1750 
r [kpc] 



Fig. 2. Deprojected (triangles) and projected (stars) X-ray gas 
temperature profiles measured in the [0.8-10.0] keV energy 
range. The data points of the projected profile are slightly shifted 
to the right to improve the clarity of the plot. The solid line shows 
the best fit function used in the total gravitational mass estima- 
tion presented in Sect. l6.1l below. 



from the maximum cluster temperature at a break radius rbr 
moving outwards and drops towards the cluster center If 
rbi is simply defined as the distance from the cluster cen- 
ter where the temperature is maximal, then rbi- - 433 ± 87 
kpc for the deprojected profile and rbi = 260 + 87 kpc 
for the projected profile, respectively. This distance corre- 
sponds to ~ 0.1 - 0.2rvir (see Sect. |6.2t , in agreement with 
works on the scaling properties of large samples of clusters of 
galaxies ("Markevitch et al."1998'; 'De Grandi & Molendil Hool 
iPiffaretti et al. 2005; Vikhlinin et al. 2005; Pratt et al. 2007]). The 
temperature decrease observed in the outer regions (~ 40% from 
Tbr to 0.5 rvh) is also consistent with the findings of these studies. 
The temperature derived from the deprojected spectral analysis 
drops from the peak value of 13.6 keV to the central minimum 
value of 9.1 keV. This is fully consistent with the typical 30% 
drop seen in temperature profiles of cool core clusters (e.g. see 
iKaastra et al.]|2004 . 

5.2. Cooling time 

The cluster RX J 1 3 4 7 is known to host a cool core 
(SchindleretiD \l99% lAIIen et all l2002bl: iGitti & SchindleJ 
2004). The centrally peaked surface brightness profile and the 
central temperature drop discussed in Sect. [3] and Sect. 15.11 re- 
spectively, are indeed signatures of the presence of a central re- 
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Table 3. Results of the deprojection analysis on annular 
MOS+pn spectra using the XSPEC pro jet model. The column 
density is fixed to the Galactic value and the normalizations are 
in units of 10"''*nenpV/47r[DA(l + z)?. The fit gives x^/dof = 
3007/2557. The data of the SE quadrant are excluded (undis- 
turbed cluster). 



Radius 
(kpc) 



0-173 
173-346 
346-519 
519-692 
692-1038 
1038-1731 



kT 
(keV) 



Q I +0.4 
'■^^-0.4 

12 6+' -^ 
13.6+^-^ 
8-6!" 

ii.r« 

7 l**-^ 



z 


norm 






(xlO-3) 


(xlO-^ cm-3) 


32+'"'' 

U.iU Q [J 

22+"-^^" 


6 02+"-"« 

9 71+0.09 
^- '-'-0.07 

1 5r''-'" 

-Q-!!9 


9T 99+0.16 
-0A6 
c q 1 +0.09 
-^•^ -0.08 

2 66+°"* 

-Q.08 




0.88!: 


1 46+°''' 

-Q.!15 


0.39!«;« 


0.81!™ 
0.54!»:«? 


69+'' ''^ 

"•"^-0.03 
"•^-'-0.02 



200 
100 
50 
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Fig. 3. Cooling time as a function of radius. 



gion where the plasma cooling time is short. In the following we 
compute the cooling time profile and the cooling radius of the 
cluster 

The cooling time is calculated as the characteristic time that 
it takes a plasma to cool isobarically through an increment of 
temperature 6T: 



^cool 



5 k6T 



2 neA(r) 



(3) 



where A(r) is the total emissivity of the plasma (the cooling 
function) and k is Boltzmann's constant. Utilizing the depro- 
jected temperature profile and the density profile from Sect. 14.31 
we can calculate the cooling time as a function of radius, which 
is shown in Fig. [3] The cooling time shows a power law behavior 
as a function of radius. We find fcooi ^ f-^-'^6±o.oi y^jjgjj ^jj g radial 
bins are used in the fit and fcooi 72±o.2i ^j^jy j-j^g 4 radial 
bins beyond 0.2 rsoo ~ 280 kpc are considered (see Sect. |6.2| be- 
low for the definition and computation of rsoo). The latter value 
agrees with recent results from the analysis in the same radial 
range of a sampl e of luminous cluste rs at z = 0.2 (IZhang et alj 
I2007h . Following Hirzan et al.1 (|2004 . we define the cooling ra- 
dius as the radius within which the gas has a cooling time less 
than 7.7 X 10'' yr, the look-back time to z = 1 for our adopted 
cosmology. With this definition, we find rcooi ~ 210 ± 10 kpc 
which corresponds to the central 36 arcsec. 

In the following analysis it is important to correct for the 
effects of the central cooling flow when measuring the char- 
acteristic temperature of the undisturbed cluster The average 
emission-weighted cluster temperature is calculated by fitting 



with a mekal model the spectrum extracted up to the outer ra- 
dius detected by our X-ray observation (5 arcmin), after excising 
the cooling region (central 35 arcsec) and the SE quadrant. We 
find a value < Tx 10.1 ± 0.7 keV. 



5.3. Entropy 

The gas entropy in groups and clusters of galaxies has recently 
received particular attention since it resulted to be a very useful 
quantity to probe the thermodynamic history of the hot gas in 
these systems. The entropy is usually defined as S - kT jn^J^, 
where T and are the deprojected electron temperature and 
density, respectively. 

In cooling core clusters the radial entropy profiles are ex- 
pected to increase mono tonically moving outvv ards, and to show 
no isentropic cores (e.g. . [McCarthy et al.l2004l) . This beh avior is 
indeed observed in nearby cooling core clusters (Piffaret fi et al.l 
2005; Pratt et al. 2006). Entropy profiles are in general well de- 
scribed by a power law. The value of the power law index scatters 
around unity, depending on the clust e r or clu ster sample used 
to derive it: for example, Ettori et al. (2002a) found 0.97 from 
Chandra data of A1795, Pratt & Arnaud (2005) derived a slope 
of 0.94 + 0.14 from scalings of th e entropy profi l es of 5 clus- 
ters observed with XMM-Newton, iPiff'aretti et all (|2005|) found 
0.95±0.02 using scaled pr ofiles of 13 cool co re clusters observed 
with XMM-Newton, and Pratt et al.' (2006) derived a slope of 
1.08 ± 0.04 (extending the sample studied in Pratt & ArnaudI 
(120051) to 10 objects). 

In Fig. |4] we show the gas entropy profile of RX J 1347 
computed from the deprojected temperature and electron den- 
sity derived in Sect. 14.31 We fit the profile with a line in log- 
log space (with errors in both coordinates) and find: log [5] = 
(1 .053 + 0.005) X log[r] + (0.01 1 + 0.010) (entropy in keV cm^ 
and radius in k pc), w hich is consistent with previous results. 
iDonahue et al.l (|2006|) recently found that the entropy profiles 
they derived from Chandra observations of 9 cool core clusters 
are better fitted by a power law plus a constant entropy pedestal 
of ^ 10 keV cm^ than by a pure power law. We performed 
similar fits and find an entropy pedestal consistent with zero. 
However, we notice that this result might be due to the lack of 
adequate spatial resolution of the entropy profile in the central 
region. 

Recent results suggest that the entropy scales with the tem- 
perature as 5 oc<Tx >''-^^, the so-called "e ntropy ramp", instead 



of the self-similar scaling S oc< Tx > dPonman et al 



2OOI 



2006h . 



iPratt & ArnaudI 120051: iPiffaretti et ail 120051: iPratt et al.1 
Here < Tx > is the mean cluster/group temperature corrected 
for the cool core effect and S is the entropy measured at some 
fraction of he virial radius (usually 0.1 x r2oo, see Sect. |6.2| be- 
low for the definition and computation of r2oo). In order to verify 
if the entropy measured in RX J 1347 follows this relation, we 
therefore adopt the scaling S oc h-'^'^{z)i< Tx > /10keV)"•''^ 
with a mean temperature for RX J 1347 equal to 10.1 ± 0.7 keV 
(see Sect. 1521). Here h^(z) = Qm(l +1^+0.^ and the factor /z-'^/^ 
comes from the scaling of the density. At 0. 1 x r2oo the scaled 
entropy is equal to 382 ± 32, 349 ± 54, and 437 + 5 1 keV cm^ for 
r2()o computed from the total mass profiles derived from model 
SO, DDgl, and NEW, respectively (see Sect. |6.1| below for the 
different models used in the total mass determination from the 
X-ray data). If instead the value Tjq™ is used (i.e., we adopt the 
size-temperature relation calibrated through numerical simula- 
tions, see Sect. 16.21 below), the normalization is 567 ± 70 keV 
cm^. The errors on these normalizations also take into account 
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Fig. 4. Entropy as a function of radius and the best fit power law 
log[5] = (1.053 ± 0.005) X log[r] + (0.011 + 0.010) (entropy in 
keV cm^ and radius in kpc). 



the uncertainty in the estimate of raoo- The normalization derived 
by adopting the size-temperature relation is in good agreement 
with the entropy normalizatio n of the 5(0.1 x r2o o)-< Tx > re- 
lation at < Tx >= lOkeV by (iPonman et al.ll2003l see their Fig. 
4). The values computed using r2oo derived from the total mass 
profiles are smaller, but still consistent within the uncertainties, 
than the values found bv lPonman et al. I (l2003l) . 

6. Mass determination 



In iGitti & Schindleii (|2004|) we presented the total mass profile 
estimated from the single yS-model. Here we perform a detailed 
study of the radial profiles of total gravitational mass and gas 
mass reconstructed by using different method s. The new values 
do not change the main conclusions in IGitti & Schindleii (12004 
but are more accurate. In this section we also present the com- 
putation of the characteristic radii quoted above. 

6.1. Total gravitational mass 

The analysis to estimate the total gravitational mass of RX J 1347 
is not limited to only one specific method, but is instead caiTied 
out by adopting different approaches. This enables us to inves- 
tigate the effects introduced by different fitting functions for the 
gas density and temperature, and different methods to derive the 
total mass from the observed gas distribution. 

The total gravitating mass distribution is calculated under the 
usual assumptions of hydrostatic equilibrium and spherical sym- 
metry by using 



Mtot(< r) = - 



kT(r) r 
Gj-mip 



d\npg{r) d\nTg{r) 



d\\\r 



d\nr 



(4) 



where G and nip are the gravitational constant and proton mass 
and ju = 0.62. A welcome property of Eq.|4]is that the total grav- 
itational mass within a sphere of radius r is determined from the 
gas density pg and temperature T„ measured at the cluster-centric 
distance r. This implies that when the gas density and tempera- 
ture are well modeled only in the radial range R\a - Rom but not 
within Rin, the mass determination is still reliable in the range 
Rin - Rout- As shown in Sect. [3] a single yS-model provides a good 
fit to the surface brightness profile in the radial range 350 kpc < 
T< 1730 kpc (model SO in Table[T]). In this case the deprojected 
gas density profile is easily computed and the total cluster mass 



is independent of the gas density central value. Since beyond 350 
kpc the temperature profile is declining, it can be well modeled 
through the poly tropic relation T cc p^^ , with 1 < y < 5/3. 
The polytropic fit to the deprojected temperature profiles gives 
in this case y - 1.23 ± 0.02 (Icr eiTor on one parameter). The 
total mass profile computed using this model is discussed be- 
low together with the results from the more sophisticated double 
y6-model. 

In order to obtain a total mass estimate for the whole ob- 
served radial range we use the double j6-model fits discussed 
in Sect. [3] (model DD and DE in Table [1]). The gas density is 
computed from the dou ble ;8-mode l surface brightness fits using 
the formulas derived in lXue & Wul ([2000i) : we assume that each 
component corresponds to a gas phase, invert Eq.|2]and compute 
the the electron number densities for the two components ne,i('") 
and the total electron number density, nf(r) using; 



1/2 



«.(0) 



njr) 
ne,i(.>') - «e,i(0) 



ne,i(r), 

2 r^'*' 



r ■ 

V c,i J 



(5) 



(6) 



(7) 



where i=l,2 and nf(0) is the central, total electron density. The 
central number densities for the two components are given by 



in which 



4;r'/2 



r(3A-l/2) 



S 0,; 
rc4 



8jl\''c,jSo.i 



1/2 



ropj) r(3A-i/2) 



r(3A) r(3/3,-i/2) 



(8) 



(9) 



where j= 1,2 and j + i. Here gi is the Gaunt factor for the compo- 
nent / and a{Ti) is the emissivity due to thermal bre msstrahlung. 
The G aunt factors are computed using the results of ISutherlandl 
(Il998h . Note that in the derivation of the equations given above 
it is assumed that each component has a constant electron tem- 
perature Tj throughout the cluster As shown in Sect. |5] the gas 
is not isothermal hence this assumption is not strictly valid. 
Nevertheless the temperature dependence of the above equation 
is fairly weak and we set rioi2 = T^ax - 13.6 keV (the maxi- 
mum of the temperature profile) and - T2oi\ - r,ni„ = 7.1 keV 
(the minimum of the temperature profile) to quantify the max- 
imum variation of the total mass estimate with temperature. 
Using the above equations and Eq.|4]we compute the mass pro- 
file for 4 cases: DDgl (model DD and Ti - Tmax, T2 - Tmin), 
DDg2 (model DD and Ti = Tmin, 7^2 - ^max), DEgl (model 
DE and T[ - T,^^^, T2 = T,^in), and DEg2 (model DE and 
Ti - Tmin, T2 = Tmax)- While the assumption of isothermal- 
ity is justified in the evaluation of the density-dependent term 
of Eq.|4]from the observed surface brightness profile, the radial 
dependence of the gas temperature must be carefully modeled, 
since the total gravitational mass varies strongly with temper- 
ature. The temperature profile in the whole observed range is 
clearly not well described by a polytropic relation and it is not 
possible to model it using a single analytical function due to the 
central temperature drop. We therefore model the profile using 
two functions joined smoothly at a cut radius 7?cut, i-e. we take 
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care that the tetnperature profile and its gradient are continuous 
across R^^i- Since the polytropic relation provides a good de- 

y— 1 

scription in the outer region, we adopt T cc p'„ as fitting func- 
tion for r > R^ai, with p computed from the double y6-model fits. 
The values obtained for the parameter y are very similar to those 
obtained when using the single j6-model. Within R^^at we choose 
to fit the temperature profile using a 5-th order polynomial with 
zero derivative at the center If the latter condition is not satis- 
fied the derived total mass density is found to be negative in the 
cluster core. We vary 7?cut and find that T^cut = 520 kpc provides 
the best model. The resulting best fit function is shown in Fig.|2l 
The total mass profiles computed from the surface brightness fits 
presented in the following are computed using this temperature 
profile modeling and will be indicated by the name of the model 
used to describe the surface brightness (see Table[T]). 

The relative difference between the mass profiles for model 
DDgl and DDg2 (DEgl and DEg2) is less than 4 % (6 %) in 
the whole observed radial range (0-1731 kpc). Models DDgl 
and DDg2, and DEg 1 and DEg2 give nearly identical results for 
r > 500 kpc. The largest difference is found between models 
DDgl and DEg2, but it is less that 15 % in the whole radial range 
and less than 5 % for r > 250. These small differences show that 
the temperature does not significantly affect the gas density de- 
termination for this massive and hot cluster, and that models DE 
an DD provide the same mass estimate for the whole radial range 
of interest. Given these results and the fact that model DD gives a 
smaller than model DE for the surface brightness modeling, 
we will discuss, in the following, only the mass profile derived 
using model DDgl. We compare the mass profiles derived from 
the double j6-model with the one from the single y6-model in the 
radial range 350-1731 kpc. In this range the relative difference 
of the mass profiles is at most 13 % (close to the innermost and 
outermost radii), but smaller than 10 % in the range 380-1500 
kpc for the four double j6-models we derived. Hence, the double 
jS-model provides estimates in good agreement with the single 
jS-model, and is of course preferred since it allows us to estimate 
the mass in the whole observed radial range, i.e. 0-1731 kpc. 
The mass profiles from the double y6-model (model DDgl) and 
the single j6-model (model SO) are plotted in Fig. |5] EiTors on 
the total gravitational masses are computed by propagating the 
Icr errors on the surface brightness and temperature profiles best 
fit parameters, and are of the order of 10 % and 20 % for the 
values derived from the single and double y6-model, respectively. 
The profile derived using the single y6-model is shown only in 
the region where it is valid, i.e. for r > 350. The mild depression 
visible around ~ 250 kpc in the mass profile derived from the 
double jS-model is due to the shape of the temperature profile in 
the inner region. 

The cluster gravitational mass can also be computed by mak- 
ing direct use of the gas temperature and gas density profiles de- 
rived from the deprojection analysis presented in Sect. 14.31 We 
invert the equation of hydrostatic equilibrium (EqU) and, using 
the three-dimensional gas density, we select the dark matter mass 
model that reproduces better the deprojected temperature pro- 
file. In the minimization the Icr errors on one single parameter 
from the spectral fits are used. For dark matter mass model, we 
consider the integrated NEW (.NavaiTO et al...l996.) dark matter 
profile: 



Mdm(< r) = 4;rr,Vc 



2oo4n(l+r/r.)-^) 



(10) 



3 ln(l -Hc)-c/(l +c) 

where pc.c = (3//?)/(87rG) is the critical density at the cluster's 
redshift. The scale radius and the concentration parameter 



c are the free parameters. The total gravitational mass within 
a sphere of radius r is given by gas plus dark matter mass 
and therefore Mtot(< r) - Mgas(< r) + Mdm(< r) in Eq.Q. 
Nevertheless, in most of the work Mtot(< r) = Mdm(< f) is used, 
i.e. the NFW profile is used to fit dark matter plus gas mass. We 
also computed the total mass profile by taking into account the 
gas mass, i.e. by adding the cumulative gas mass profile to the 
best-fitting NFW profile, and found little difference between the 
two profiles. The best-fit parameters are rs = 722 + 112 kpc and 
c - 3.20 + 0.30 (eiTors are RMS of the Icr joint confidence lim- 
its), with;i'^j|^ = 6.7 for 4 degrees of freedom. Our best-fit NFW 
profile is shown in Fig.|5] From the set of (c, parameters ac- 
ceptable at 1 cr we compute, for each radius, the maximum and 
minimum value of the total mass and hence its upper and lower 
errors. These are of the order of 10 %. From a visual inspection 
of Fig. |5] one can note that the NFW mass profile is lower than 
the double p estimate for r < 1 150 kpc and higher at larger radii. 
The discrepancy within 1 150 kpc is due to the fact that our best- 
fit NFW profile tends to underestimate the temperature in this 
range. The relative difference between the NFW and the double 
/? mass profiles is -38 % (underestimate) at r = 500 but decreas- 
ing towards the center, and increases almost linearly to 30 % 
(overestimate) at r = 1731 kpc. The fairly low concentration 
parameter c, compared to the p redictions of numerical simula- 
tions (e.g, iMaccio' et al.ll200"6l) . and the goodness of our NFW 
fit might indicate that our temperature profile is not enough spa- 
tially resolved in the central region of the cluster for this kind of 
mass determination method. While the mass determination from 
the double yS-model may therefore be preferred, we present val- 
ues also from the NFW fitting for completeness. 

6.2. Virial radius and scaling reiations 

In this section we determine the characteristic radii ta used in 
Sects. Island 15.31 For the various mass profiles we compute r^, 
the radius within which the mean interior density is A times the 
critical value, by using 



A = 



3Mtot(< ta) 



(11) 



For the cosmology adopted here the virial radius is given by 
'"vii - r^, with A = 178-1- 82;ic - 39x^ and where x - Q(z) - 1 and 
Q(z) = 0.3 (1 + z)V(0.3 (l+z?+ 0.7) (Bryan & Norman 1998i). 
Thus for RX J1347 A = 135. We also compute Mtot(< /"a) and 
Mgas(< ta) for various overdensities: A = 2500, 1000,500,200. 
The results obtained from the overdensity profiles calculated 
from the double yS-model (DDgl) and NFW fit are reported in 
Tables 

The size-temperature relation ta oc y< Tx > predicted by 
self-similarity allows an estimate of ta from the mean clus- 
ter temperature alone, provided that its normalization is known 
from numerical simulations. We compute the normalization for 
the cosmolog y adopted here by interpolating the values given in 
lEvrard et al.l (11996). For the mean cluster temperature < Tx >- 



0.7 keV we derive the characteristic radii r^™, finding 

Sim _ OQfZ , ir\ „„j ^Sim 



10.1 



'2500 = ± 30 kpc and r^™ = 3197 + 107 kpc. From our 
X-ray analysis we find r25oo (734 ± 34,729 + 63,608 ± 20) 
kpc and r^i, = (2378 + 76,2241 + 189,2639 + 108) kpc when 
using in Eq.[TT]the mass profile derived from model (SO, DDgl, 
NFW), respectively. These values are consistent with the size- 
tempera ture relation derived from observations of nearby relaxed 
clusters dArnaud et al.ll2005l) . By comparing the above values we 
note that the estimates from the X-ray analysis are systemati- 
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Table 4. Characteristic radii ta, total mass Mtot and gas mass Mgas for various overdensities A derived from the double yS-model 
(DDgl) and NFW fits (Icr errors in parentheses). The masses are estimated within ta. As discussed in Sect. 16.11 results from the 
double /3-model are generally more reliable. 



A 


'A.DDsl 


A^tot.DDgl 


M.,as.DDgl 


A.NFW 


Mol,NFW 






(kpc) 


(10"*Mo) 


(l6"*Mo) 


(kpc) 


(10"*Mo) 


(10»Mo) 


200 


1957.2 (183.2) 


11.00 (2.78) 


3.34 (0.11) 


2286.7 (110.8) 


17.86 (2.07) 


3.82 (0.18) 


500 


1387.2 (123.9) 


9.77 (2.30) 


2.39 (0.05) 


1479.2 (71.23) 


11.85 (1.13) 


2.55 (0.06) 


1000 


1063.3 (91.00) 


8.80(1.91) 


1.79 (0.03) 


1029.3 (43.08) 


7.99 (0.62) 


1.73 (0.03) 


2500 


729.3 (63.2) 


7.10(1.40) 


1.15 (0.02) 


608.1 (19.8) 


4.12 (0.22) 


0.91 (0.01) 
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Fig. 5. Integrated three dimensional total mass profiles, with er- 
rors, derived from the double /3-model (model DDgl, solid), 
single y6-model (model SO, dashed), and the NFW model (dot- 
dashed). The dotted line shows the cumulative gas mass profile. 
See text for details. 



cally lower than the ones predicted from the size-temperature 
relation calibrated by means of numerical simulations. It is not 
surprising that we find a smaller discrepancy for r250() than ryu, 
as its determination does not require extrapolation of the ob- 
served mass profile. Thi s is in agreement with results for other 
individual clusters (e.g., Gitti et al. 2007) and studies of clus- 
ter samples (Sanderson et al. 2003; Piffaretti et al. 2005). The 
largest discrepancy is found for ryir and in poor, cool clusters. In 
these systems the impact of additional, non-gravitational heating 
is most pronounced, as the extra energy required to account for 
their observed properties is comparable to their thermal energy 
(Ponman et al. 1996; Tozzi & Norman 2001). The observed dis- 
crepancy is also related to the cluster total mass determination. 
In this context it is interesting to note that recent results from 



numerical simulations indicate that the total mass of simulated 
clusters estimated through the X-ray approach is lower that the 
true one due to gas bulk motions (i.e. deviation from the hydro- 
static equihbrium ) and the coraplex thermal structure of the gas 
dRasia et alJl2006tlNagai etalllioOTh . 

A self-similar scaling relation between Mtot and < Tx > at 
a given overdensity is predicted in the form Mtot Tx >^^^. 
Various observational studies have found different and some- 
time conflicting results regarding the slope and normalization of 
the M-T rel ation (e . g.,[A llen et al. 2001; Finoguenov et al. 200ll 
lEttori et all l2002bl: [Sanderson et alj i2003i; |Amaud et al. 20051 
and re ferences therein). The relation derived by Arnaud e t al] 
(I2005h for a sub-sample of six relaxed clusters hotter than 3.5 
keV observed with XMM-Newton is consistent with the standard 
self-similar expectation, following the relation: 

J ^> 1.51±0.11 

Kz)M2500 = (1.79 ± 0.06) X lO'^^Mo (^^) (12) 

This result is i n agreement with Chandra observations 
(lAllen et al.ll200lh . In the case of RX J1347, Eq.[l2] turns into 
an estimate of M2500 - (4.07 + 0.46) x IO'^^Mq. By consider- 
ing the whole XMM-Newton sample (ten clusters in the temper- 
ature range [2-9] ke V), the relation steepens with a slope ~ 1 .70 
(lArnaud et al.ll2005l) indicating a breaking of self-similarity. In 
this case we estimate M25oo,DDgi = (4.39 ± 0.35) x IO'^^Mq. The 
mass estimate that we derive at the overdensity A - 2500 differs 
strongly depending on the model adopted (see Table |4|. From 
model DDgl we estimate Masoo.oDgi = (7.10 ± 1.40) x IO'^^Mq, 
which is much higher than the prediction of the M-T relation. 
The mass estimate of M25()o,nfw = (4.12 ± 0.22) x IO'^^Mq as de- 
rived from the best-fitting NFW profile is instead in good agree- 
ment with the M-T relation, although the large error bars prevent 
us from distinguishing between a self-similar or steeper relation. 

6.3. Gas mass and gas mass fraction 

From the results of the deprojected spectral analysis we compute 
the cumulative gas mass profile Mgas(< r), thus obtaining values 
for the 6 bins used in in the spectral analysis. In order to derive 
better estimates when an extrapolation of the gas mass beyond 
Rout is needed, we compute the gas mass profile using the radial 
gas density profile derived from the best fit parameters of the 
double j6-model (model DDgl) of the surface brightness profile. 
The normalization of the latter is fixed using the gas density pro- 
file from the spectral analysis. The resulting gas mass profile is 
shown in Fig.|5] When Mga,s(< r) is evaluated within 7?out we use 
the binned profile and spline interpolation, which in this radial 
range provides values consistent with the ones computed using 
the results from the double j6-model. 

The gas mass fraction /gas is defined as the ratio of the total 
gas mass to the total gravitating mass within a fixed volume. We 
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measure /gas,250() — 0.162 ± 0.036 from the mass profiles derived 
from the double y6-model fit (model DDgl). This value is close to 
the global baryon fraction in the Universe, cons trained by CMB 
observations to be Q^/Q m = 0.175 ± 0.023 dReadhead et alj 
120041; ISpergel et al.ll2003h . and is higher than the average value 
deriv ed in a number of p revious measurement s with Chandra 
(e.g., lAUen et al.1 l2002al: IViMilinin etall l2006 i However, we 
note that a general trend of increasing /ga s with cluster tem pera- 
ture (hence mass) has been observed ( Vikh linin et al.ll2006h . The 
high central gas mass fraction measured here is consistent with 
this tendency, as RX J 1347 is a hot, massive cluster. 



7. Comparison with previous work 

In this section we compare the most relevant total and gas mass 
estimates for RX J 1347 found in literature with our results. The 
values in literature are converted to the cosmology adopted here 
before the comparison. 



7. 1 . Comparison with X-ray studies 

Using combi n ed R OSAT and ASCA observations 
ISchindler et al] (1 19971) derived M,ot = 1.11 x 1O"*M0, 
Mft = 4.93 X IQ^'^Mq, and Mft = 1.45 x lO'^^M© within 204, 
850, 2550 kpc, respectively. These values were derived assum- 
ing isothermality and the eiTor coming from the uncertainty 
on the global temperature is of the order of 10%-15%, as we 
estimated fro m the plot showing the profile of the integrated 
total mass (see ISchindler et an [T997. Fig. 6). We find, for model 
(SO, DDgl), Mtot = (1.14 + 0.14,0.93 + 0.17) x lO'-^M©, 
Mtot = (8.10 + 0.84,7.85 ± 1.60) x IO'^^Mq, and 
M,ot = (1.47 + 0.15,1.22 + 0.32) x lO'^M© within 204, 
850, 2550 kpc, respectively. While we find a significant mis- 
match at 850 kpc, the results are in reasonably good agreement 
at small and large radii, in particular considering the errors 
and the different assumptions adopted in the mass d e termi- 
nation. For the cumulative gas mass ISchindler et al.l (Il997h 
found Mgas = 1.33 X lO^^M© and Mgas = 5.93 X lO^^^M© 
within 850, 2550 kpc respectively, while our values are 
Mgas = (1.39±0.02)x10'4Mo and Mgas = (4.32+0.17)xlO''*Mo 
within 850, 2550 kpc, respectively. While the v alues at 850 kpc 
are co nsistent, the large value found at 2550 by ISchindler et alj 
(199% is very likely due to the narrower radial range probed 
by their observation. As shown in Sect. [3] (see Table [1]), the gas 
density steepens in the outer region. As a result, the gas mass 
derived from a single y6-model fit to a narrow central region 
and extrapolated to large radii is biased high. In comparing 
the results, we should also bear in mind that the analysis 
presented by ISchindler et al.l (119971) is performed on the full 
360° data, as the hot enhancement in the SB quadrant has been 
discovered o nly subsequently with Chandra and XMM-N ewton 
observations dAllen et al.l l2002b'; 'Gitt i & Schindleiil2004 . 

We c ompare our best-fitting NFW profile with the one de- 
rived bv lAUen et alj (l2002b) from Chandra data. The two pro- 
files are consistent, with a relative difference ranging from 15% 
to 30% depending on the radial range considered. As a gen- 
eral trend, our profile results lower in the inner region (in- 
side ~600 kpc) and higher i n the outer r egion ( outside ~1000 
kpc) than the one d erived bv lAllen et al.l (l2002bl) . In particular, 
lAllen et aTl (l2002b l) find an integrated mass within the virial ra- 
dius of their best-fitting NFW mass profile of Mtot(< 2Mpc) - 
(1.95^q 7q) X 1O'^M0, which is in fairly good agreement with the 
value that we measure: M,ot(< 2Mpc) - (1.59+[5 ]g) x lO'^M©. 



When considering the mass profile derived from the double j8- 
model (DDgl), which at large distances is lower than the one 
derived from the NFW fit (see Fig.|5]), we find Mtot(< 2Mpc) = 
(1.1 1+15 27) ^ l O^^MfD. Thi s value is fairly low compared to the 
value found bv lAllen et al.l (l2002bl) . but still consistent consider- 
ing the errors on the mass estimates at this large distance. 

lEttori et al] (|2004|) derive from Chandra data estimates of 
Mtot = (8.94+0.80)xl0i'*Mo andMgas = (1.81±O.O8)xlO"*M0 
within 1368 kpc, which corresponds to r^oo in their work. While 
our value Mtot(< 1368kpc) = (10.84 ± 1.11,9.72 + 2.27) x 
lO'^^Mo (for model SO and DDgl, respectively) agrees with 
the Chandra estimate, we find a larger value for the gas mass: 
Mgas(< 1368kpc) = (2.35 + 0.05) x 10''*Mo. The discrepancy 
might be related to the different approaches adopted for the cal- 
culation. We estimate the gas mass directly from the density pro- 
file derived from the deproj ected spec t ral an alysis (Sect. 14. 3t . 
The gas mass computed by Ettori et al.l ( |2004|) is derived by es- 
timating the central electron density from the combination of 
the best-fit results of the spectral and imaging analyses (namely 
the normalization of the thermal spectrum and the parameters 
of the singl e ,6-mo del). In particular, the low value measured by 
lEttori et all (|2004|) might be biased by an underestimate of the 
central density due to a possible undersampling of the cluster 
luminosity within the radius where the thermal spectrum is ex- 
tracted, which coiTesponds to only ~ 0.4 times the radius where 
the X-ray emission is detectable in the Chandra data. 



7.2. Comparison witli dynamical estimates 

Using the virial approach, ICohen & KneibI (l2002h derive Mtot = 
(6.92^j-^g) X IO^'^Mq within 597 kpc from galaxies velocity dis- 
persion measurements. Within this radius we consistently find 
Mtot = (6.06+0.61, 6.05+1.14)xlO''*Moformodel (SO, DDgl). 



7.3. Comparison witli gravitational lensing 

Since the gravitational lensing analysis measures the projected 
total mass distribution, in order to compare consistently the re- 
sults from the X-ray and lensing techniques we project along the 
line of sight the cumulative 3D mass profiles Mtot(< r) derived 
in Sect.m thus obtaining M^^°'(< r). 

From a weak lensing investigation, iFischer & TvsonI (|1997[) 
derive Mf^ = (9.35±2.55)x1O14M0 within 850 kpc. We find, in 

good agreement, Mf^ = (9.95 ± 1.03,9.17 + 2.13.7 9) x lO'^M© 
for model (SO, DDgl). Note that lFischer & TvsonI fl 997) com- 
pare their mass measurement Mjl°^ wit hin 850 kpc with the Mtot 
value quoted in ISchindler et alj d 19971) a nd find a large discrep- 
ancy. As pointed out by Sahu et alT i 19981) . the two mass determi- 
nations are in agreement if the correct quantities are compared. 

In the strong lensing analysis bv' Sahu et alj dl99 8), M\^°^ - 
5.36 X 10'"* M0 is measured within 204 kpc, the cluster-centric 
distance of the arcs. Within this projected distance we find 
Mf^j = (2.38 ± 0.27,2.05 + 0.37) x lO'^^M© for model (SO, 
DDgl). Although this discrepancy might be due to the fact that 
we excise the perturbed region of the cluster, we note that such 
a large mismatch between the masses dete rmined from X- rays 
and strong lensing is commonly found (see lWu et alJll998l and 
references therein). In the inner core of clusters, where strong 
lensing occurs, the physics of the ICM may be complicated by 
the interaction with the central AGN. The central cluster region 
is thus poorly described by the usual simple models used in the 
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Fig. 6. Ratio of the lensing to the (projected) X-ray mass profile 
for different X-ray mass estimates. The line styles are the same 
as in Fig.|5] The reported errors are those coming from the X-ray 
mass determinations. 

X-ray methods which rely on the assumptions of spherical sym- 
metry and hydrostatic equilibrium. 

We com pare our total mass determination wi th the lens- 
ing resu lts of iBradac et a l.l (2005b). The results of Bradac et al.J 
( l2005bl) are obtained using a mass reconstruction method which 
combines strong and weak gravitational lensing data and effec- 
tively breaks the mass-sheet degeneracy (Bradac et al. 2005a|. In 
Fig. |6]we show the X-ray to lensing mass ratio Miensing /A^x-ray 
as a function of radius up to ~ 670 kpc, the limiting radius of 
the lensing study. From a visual inspection of this figure it is 
clear that there is lack of agreement between the X-ray and lens- 
ing mass estimates. Only in the central region the X-ray mass is 
marginally consistent with the lensing mass. The mass ratios in- 
crease with radius and tend to approach a constant value at large 
radii. At 600 kpc the ratio is 2.07, 2.17, and 2.45 for the X-ray 
mass estimated using the single /3-model (SO), double y6-model 
(DDgl), and NFW model, respectively. We stress that the same 
discrepancy is found when we compare our mass profiles with 
a corrected mass profile computed from the lensing map where 
the SB quadrant, which contains the hot X-ray subclump, is ex- 
cluded. 

7.4. Comparison with tlie Sunyaev-Zei'dovich (SZ) effect 

Through the SZ effect, iPointecouteau et aTl ( l200ll) measure the 
gas mass of R X J1347. They c ompa re their results with the X- 
ray results of ISchindler et al.l (1 19971) . finding good agreement. 
Within 74 arcsec= 427 kpc the SZ estimate is Mg^s - (4.7 + 
0.4) X 1O^^M0 in agreement with our value, Mgas - (5.5 + 0.1) x 
IO'^Mq. 

8. The cool core 

As shown in Sect. |5] there is no evidence for very low temper- 
ature gas in the core of RX J 1347, suggesting that the descrip- 
tion of the inner region of this cluster by means of a standard 
cooling flow model is not appropriate. The spectral analysis in 
iGitti & Schindled(l2004. ) shows that if the cool core in RX J1347 
is fitted with an empirical cooling flow model where the lowest 
temperature is left as a free parameter, very tight constraints on 
the existence of a minimum temperature (~ 2 keV) are found. 
This situation is common for cool core clusters and it has be- 
come clear that the gas with short cooling time at the center of 



these objects must be prevented from cooling below the observed 
central temperature minimum. The most appealing mechanism 
is heating by AGN because it is strongly motivated by observa- 
tions. Central AGNs with strong r adio activity are found in the 
majority of cool core clusters (e.g. lBurrisll 19901; [Ball et al.ll 19931) 
and powerful i nteractions of the radio sources with the ICM are 
observed (e.g.. lBirzan et alj|2004t iRaffertv et al.ll200 6', and ref- 
erences therein). The presence of a central AGN in RX J 1347 is 
indicated by the NRAO VLA Sky Survey (NVSS), that shows a 
strong central source along with some hint of a possible extended 
emission. However, the resolution and sensitivity of the NVSS 
are not sufficient to study the characteristics of the central source 
and establish the existence of diffuse emission. We obtained new 
VLA data in order to further inve stigate the nature and properties 
of the radio source in RX J1347 ( iGitti et al.ll2007 l). 

In order to explore the heating by AGN, we adopt the model 
developed by iRuszkowski & BegelmanI (l2002l hereafter effer- 
vescent heating). The details of the model and the procedure 
adopted to estimate the AGN parameters fro m the observed tem- 
perature and density profiles are given in iPiffaretti & Kaastral 
(2006). Here we simply summarize the essential elements. In 
the effervescent heating scenario the central AGN is assumed 
to inject buoyant bubbles into the ICM, which heat the ambi- 
ent medium by doing PdV work as they rise and expand adia- 
batically. In addition, besides being essential in stabilizing the 
model, thermal conduction transports energy from the hotter, 
outer region to the central region. Unfortunately its efficiency 
is poorly known since it depends on magnetic fields and mod- 
els with different fractions /c of the Spitzer rate are studied. For 
a fixed /c between and 1/3 (the maximum for magnetized a 
plasma) the contribution of heat conduction as a function of ra- 
dius is known since the temperature gradient is estimated from 
the deprojected temperature profile. We note that if one assumes 
that heat conduction alone balances radiative losses, then its ef- 
ficiency would be much larger that 1/3 of the Spitzer rate and 
therefore unrealistic. The raising entropy profile (in Sect. 15.31 ) 
indicates that convection is not operating on the scales that we 
are able to resolve and is therefore not included in the model. The 
extra heating profile resulting from subtracting the heat conduc- 
tion yield from the ICM emissivity is then assumed to be bal- 
anced by the AGN heating function: 

where p is the ICM pressure (po some reference value) and jh 
the adiabatic index of the buoyant bubbles, which is fixed to 
4/3 (i.e., relativistic bubbles). Fitting Eq.(fT3b to the extra heat- 
ing profile provides the AGN parameters L (the time-averaged 
luminosity) and (the scale radius where the bubbles start ris- 
ing in the ICM). Only if 0.10 < fc < 0.27 the fitting pro- 
vides meaningful results. For fc - 0.27 the AGN parameters are 
L = 7.45 X 10^^ erg s"' and ro = 4 kpc. As we decrease f both 
AGN parameters increase monotonically and reach the maxi- 
mum at fc - 0.10, for which we find L - 10.11 x 10"*^ ergs"' 
and ro = 29 kpc. The trend of the AGN parameters with fc indi- 
cates that, in the framework of the effervescent heating scenario, 
heat conduction and AGN heating cooperate in quenching radia- 
tive cooling. The inferred AGN time-averaged luminosity lies 
therefore in a quite small range (7.45 - 10.1 1 x 10^^ erg s"'), and 
is larger but comparable to the cluster luminosity in the energy 
range [2.0-10.0] keV (Lx = 6.2 ± 0.2 x 10'*^ erg s"'). The model 
with fc = 0.22 is the one with the smallest reduced and in this 
case L = 8.32 x lO'*^ erg s"' and ro = 13 kpc. 
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The effervescent heating model applied to RX J 1 347 predicts 
that the scale where the bubbles start rising in the ICM is in the 
range 4-29 kpc. The observed extens ion of the AGN jets shoul d 
be of the same order of magnitude (iPiffaretti & Kaastrall2006l) . 
Interestingly, the first results from 1.4 GHz VLA observations 
of the central region of RX J 1347 show hints of faint structures 
emanating from the discrete radio source out to ~ 20 kpc from 
the center (Gitti et al. 2007). A comparison between the derived 
luminosity L with the observed AGN luminosity is unfortunately 
not possible. In fact, in the framework of the effervescent heat- 
ing model, the derived AGN luminosity is a time-averaged total 
AGN power and a fair comparison is possible only if the total 
jet power is estimated (X-ray and radio powers are known to be 
poor tracers of the total AGN po wer). At present, th is was done 
only for M87 in the Virgo cluster dOwen et al.ll2000l) . 



9. Summary 

As indicated by previous studies ([Allen et al.l l2002bl: 



iPointecouteau et al.ll200lt iGitti & Schindleil l2004l)T the cluster 
RX J 1347 shows both the signatures of strong cooling flow and 
subcluster merger, that are rarely observed in the same system. 

We analyze the data excluding the SE quadrant, where the 
presence of a hot X-ray subclump is suggesting that a minor 
merger has recently occurred or is still going on, and find that: 

- the features (shape, normalization, scaling properties) of 
density, temperature, entropy, and cooling time profiles are 
fully consistent with those of relaxed, cool core clusters, with 
no indications of perturbations that may originate from the 
disturbed region of the cluster. 

The usual assumptions of hydrostatic equilibrium and spher- 
ical symmetry can therefore be adopted when analyzing the data 
with the SE quadrant excluded. This allows us to perform a de- 
tailed mass reconstruction by starting from the temperature and 
density profiles derived from the X-ray analysis. We find that: 

- the total mass profiles computed from a double and single 
jS-model for the surface brightness give consistent results if 
the cool core is excised in the latter case; 

- there is a reasonably good agreement between the total mass 
profile estimated from a double y6-model and from the as- 
sumption of a NEW profile. The differences between these 
estimates might come from a poor spatial resolution of the 
density and temperature profiles in the central region, which 
could bias the NEW method; 

- the characteristic radii ta computed from the mass profile 
are in agreement with the observed size-temperature rela- 
tion, although they are systematically lower than those de- 
rived by calibrating the relation with numerical simulations. 
The mass estimated from the NEW profile is in agreement 
with the observed mass-temperature relation, whereas that 
derived from the double j6-model profile is a factor ~ 1.7 
higher. 

We compare our gas and total mass estimates with previous 
work and find that: 

- our estimates of gas and total mass are generally in good 
agreement with those from previous X-ray, dynamical, weak 
lensing and SZ studies; 

- a discrepancy of a factor ~ 2 between strong lensing and 
X-ray mass determinations is confirmed; 



- there is a large discrepancy at all radii between our total mass 
estimate and the mass reconstructed through the combination 
of both strong and weak lensing. 

We study the AGN heating in RX J 1347 by applying the 
effervescent heating model. We find support to the picture that 
AGN heating and heat conduction cooperate in balancing radia- 
tive losses. Our predictions concerning the extension of the AGN 
jets in RX J 1347 are consistent with recent radio observations of 
the radio source at the cluster center 
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